Bright and Dark Solitons in a Periodically Attractive and Expulsive Potential 
with Nonlinearities Modulated in Space and Time 
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This work deals with soUton solutions of the nonlinear Schroedinger equation with a diversity of nonlinear- 
ities. We solve the equation in a potential which oscillates in time between attractive and expulsive behavior, 
in the presence of nonlinearities which are modulated in space and time. Despite the presence of the periodi- 
cally expulsive behavior of the potential, the results show that the nonlinear equation can support a diversity of 
localized excitations of the bright and dark type which are remarkably robust. 
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Solitons are known to appear in the nonlinear Schroedinger 
equation (NLSE) in a diversity of situations. Usually, the non- 
linear interactions are of cubic nature, but there are systems 
which engender cubic, quintic, cubic plus quintic, and 
other forms of nonlinearities of current interest to the study of 
Bose-Einstein condensates (BECs). The presence of BEC of 
trapped atoms was experimentally realized in S 0] . Not 
too much later, dark solitons in BEC were formed with repul- 
Rb atoms in and bright solitons were generated 



sive 
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for attractive ^Li atoms in iSllOi]. Along the years, many in- 
teresting advancements and motivations have appeared, as we 
can see, for instance, in Refs. lfll[[l2ll . 

There are several distinct lines of investigation of the 
NLSE. An interesting possibility was initiated in [13], with 
the introduction of a procedure to deal with the NLSE with 
varying coefficients. It is based on a similarity transformation, 
which transforms nonautonomous NLSE into stationary equa- 
tion which is easier to solve 0, . In particular, in the re- 
cent work [15] one deals with the cubic NLSE where both the 
trapping potential and the cubic nonlinearity may now vary, 
being functions of both space and time. This procedure was 
recently extended to the case of cubic and quintic nonlineari- 
ties in 0]. 

Another line of investigation of the NLSE was initiated in 
[\J], in which one considers the possibility of changing the 
potential from the usual attractive harmonic form to the well 
distinct expulsive behavior, tuned through the Feshbach reso- 
nance mechanism. This possibility has initiated several inves- 
tigations with expulsive potentials, with a variety of motiva- 
tions of current interest ||18[|T9[|2q[|2T|]. 



Since it is believed that atomic matter nonlinear excita- 
tions are of importance for the development of applications 
of BECs, we think it is of interest to develop new procedures 
which allow us to investigate the presence of bright and dark 
solitons in realistic models. With this motivation, in this Letter 
we shall join together the two possibilities above, that is, we 
shall consider the NLSE with space- and time-dependent non- 
linearities in the presence of space-dependent potential which 
oscillates periodically in time from attractive to expulsive be- 
havior. The main results are that despite the expulsive behav- 
ior of the potential, the system is still able to support soliton 



excitations of both the bright and dark type, with a diversity 
of specific features, driven by the parameters which specify 
the system. The results are robust since the solutions run in 
time with periodic or quasiperiodic features without loosing 
the localized behavior in space. 

We shall implement the procedure focusing attention on the 
NLSE of the general form 



(1) 



where v{x^t) is the potential and P(|^p) is the polynomial 
function 
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The coefficients g2n-\-i{xjt) describe the nonlinearities 
present in the system. Here we are using standard notation, 
in which x, t and all other quantities are dimensionless. 

We assume validity of the one dimensional approach, which 
is an approximation of good practical use if one works under 
the restrictions <C <C A <C a|| , for ag being the scatter- 
ing length, A the spatial scale of the wave packet, and a± and 
a|| being the characteristics transverse and longitudinal trap 
lengths, respectively. We also remark that issues concerning 
parametric amplification of elementary excitations due to pe- 
riodic modulation of the trapping potential and nonlinearities 
were studied before in [22] . Moreover, the potential and non- 
linearities to be used in this work are typical of BECs, so that 
the results obtained below may stimulate new experiments in 
the field. 

The main idea of the procedure is to transform the nonau- 
tonomous NLSE ([T]) into the stationary equation 



(3) 



where fi is the eigenvalue of the nonUnear equation above, and 
g(|$|2) has the form 
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G2„+l|*|'", 
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where the coefficients G2n+i are now real and constant pa 
rameters. To do this, we write the solution of ([B as llsilTsI] 



(5) 



The substitution of (O into ([B leads to ©, but now we have 
to have 



dp d{p^{d(t)/dx)) 



dx 



= 0, 



dt dx dx ' dx 



(6a) 
(6b) 



Next, we introduce a new function i{x^t) such that 
C(x, t) = F(^(x, t)). If we make the simple choice ^(x, t) = 
^/x(0' we can determine the width of the localized solution 
as x(^)- With this, we obtain the new equations 



524+1 (a;, = 



2i+l 



1,2,3, 



(7a) 
(7b) 

(7c) 
(7d) 



The choice of F(^) is to be done consistently. 

To illustrate the procedure with examples of interest, let 
us now focus attention to the case of specific nonlinearities. 
There are many possibilities to choose the nonlinearities, and 
here we suppose that the cubic nonlinearity is given explicitly 
by 



^3 



(x, t) = sgn(G3) X"' [1 + Aexp(-^2)]^ 



(8) 



where A is real parameter which controls the behavior of the 
cubic nonlinearity - here we will use A = 0.5 where required. 
With this choice we have 

p{x,t) = abs(G3)^/^x-'/'[l + Aexp(-e2)]-i/2, (9) 

and the function F(^) has to be obtained from dF/d^ = 
abs(G3)"^/^[l + Aexp(-^2)]. The potential in Eq. ^ has 
to have the form 



v{x,t) = cj(t)x^ + /(x, t)— /i 



[1 + Aexp(-^ 
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abs(G3)2/3x2 
where uo is time-dependent function such that 



(10) 



(11) 



and f{x^t) is an awkward expression which comes from the 
first term in (TTcl) . The other nonlinearities are given according 
to Eq. (iTdl) . We remark that the last two terms in the potential 
JTOl) contribute most significantly at small distances, as we 
show explicitly below. 
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FIG. 1: (Color on line) Plots of io(t) in black (dashed) and blue 
(solid) lines, for a = 0.1 and /3 = 0, and for a = P = 0.1, re- 
spectively. 



There are several possibilities to choose the width of the 
localized excitation. In the present work we consider 



1 



1 + [1 + a sin(t) + (3 sm{V2t)]^ 



(12) 



where a and f3 are real parameters. This choice leads to u;{t) 
according to Eq. ([TT1) and in Fig. 1 we plot u;{t) for a = 0.1 
and (3 = 0, and for a = (3 = 0.1, to illustrate its periodic and 
quasiperiodic features. 

The above results are general and can be used to investigate 
systems of current interest. Let us first consider the NLSE 
in the case of cubic nonlinearity alone. This is the case of a 
BEC where the two-body scattering is strong enough to al- 
low neglecting the other higher order effects. For instance, in 
a BEC system of atomic species with naturally large interac- 
tions such as ^^Rb [4] and ^^Na [6] have properties favorable 
to be described by the cubic nonlinearity above. This leads to 
the stationary equation 



G3|$(C)|'*(C)- (13) 



We now explore the presence of bright and dark solitons. 
We first take p= — 1 e G3 = —2 to get to the potential 

t)=u;(t)x2+/(a;, t) +2-'/'x"' [1 + A exp(-C')]2, (14) 

with g3{x,t) given by In Fig. 2(a) we plot the potential 
(O to show how it behaves in space and time. Moreover, for 
/X = — 1 e G3 = —2 the equation (fT3l i has the solution 



^(C)=sech(C), 



(15) 



which is of the bright type. The wave function which solves 
the corresponding NLSE gets to the form 



^(x,t) = 2i/V^sech(C){x[l + Aexp(-^2)]} 



-1/2 



(16) 



where (j) = (j){x^t) is real, obtained via the equation (TTbl) . 

To get to the case of dark solitons, we consider /i = 2 and 
G3 = 2. In this case the potential becomes 

^(x,t)=c^(t)xV(^,^)-2^/V[l + Aexp(-e2)]^(17) 
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FIG. 2: (Color on line) Plots of t) in the presence of cubic non- 
linearity, with A = 0.5 and for (jj{t) periodic, with a = 0.1 and 
/3 = 0, as shown in Fig. 1. The plot (a) is of ([T4]) for bright solitons, 
and the plot (b) is of JTtI ) for dark solitons. 
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FIG. 3: (Color on line) Plots of v{x,t) in the presence of cubic 
nonlinearity, as given in the former Fig. 2, but now in the range 
— 2 < X < 2, to expose the behavior of the potential at small dis- 
tances. 

In Fig. 2(b), we plot the potential (fTTl) to illustrate how it be- 
haves in space and time. If we use ji = 2 and G3 = 2 in 
Eq. ([T3I) , we get the solution 

$(C)=tanh(C), (18) 

which is of the dark type. The wave function which solves the 
corresponding NLSE gets to the form 

= 2i/Vnanh(C){x[l + Aexp(-e2)]}-i/2, (19) 

with (j) = (j){x, t) real, obtained via the equation dTbl) . 

We note from Fig. 2 that the potentials oscillate from at- 
tractive to expulsive behavior, and that they show some struc- 
ture at small distances. To better see this, in Fig. 3 we plot 
the same potentials, but now we look at the behavior at small 
distances. We see that there are attractive and expulsive struc- 
tures at small distances, but they do not contribute importantly 
to the soliton formation. This conclusion follows from the fact 
that if we change A —A in ([8]), the structure at small dis- 
tances changes significantly, but the soliton solutions maintain 
the same form, with changes mostly in the amplitude. The 
other potentials show similar behavior, and the issue will be 
further explored in another work [23]. 

We now focus attention to the bright and dark solitons in 
the presence of cubic nonlinearity. In Figs. 4(a) and 4(b) we 
depict according to ([T6l) , for uj{t) periodic and quasiperi- 
odic, respectively. In Figs. 5(a) and 5(b) we depict ac- 
cording to ([T9l) , for u;{t) periodic and quasiperiodic, respec- 
tively. We note that the solutions evolve in time periodically 
or quasiperiodically, depending on the way we choose x{t) in 
([T2I) . with no further changes, showing that they are stable lo- 
calized excitations which we name periodic and quasiperiodic 
bright and dark solitons of the corresponding cubic NLSE. 
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FIG. 4: (Color on line) Plots of |^(x,t)|2 with A = 0.5 for the 
periodic (a) or quasiperiodic (b) bright soliton, with a = 0.1 and f3 = 
Oora = /3 = 0.1, respectively, in the case of cubic nonlinearity. 
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FIG. 5: (Color on line) Plots of |^(x,t)|^ with A = 0.5 for the 
periodic (a) or quasiperiodic (b) dark soliton, with a = 0.1 and f3 — 
Oora = /3 = 0.1, respectively, in the case of cubic nonlinearity. 

We now consider the case of cubic and quintic nonlineari- 
ties. This case is also of interest, since the presence of quintic 
nonlinearity is in general due to three-body interactions, and 
the cubic and quintic nonlinearities may be used to map the 
case where the two-body scattering is somehow weakened in 
a way such that the three-body effects cannot be neglected 
anymore [24] . Moreover, it may be also used to help control 
deviation of the trapped condensate from one-dimensionality 
1I25I] . Here we use the same cubic nonlinearity ([8]), but now 
the quintic nonlinearity has the form 

55 = G5 abs(G3)-"/'[l + Aexp(-e')]^ (20) 

In this case, the function obeys 

= + G,m)?m + ^51^(0^^(0. 

(21) 

The investigation is similar to the case considered in [16]. We 
first consider /i = 0, and take the cubic and quintic parameters 
as G3 = 2 and G5 = —3. In this case we get the potential 

v{x, t) = uo{t)x^ + /(x, t), (22) 

with g^{x^ t) and ^5(x, t) given by dH) and ([20b , respectively. 
Moreover, for /i = 0, and for G3 = 2 and G5 = —3, we can 
write Eq. © in the form 

— -^ = 2^3(0-3^^(0. (23) 



The solution is given by ^{Q = 1/ yl + C^. It has the bell 
shape form, and it is of the bright type. In this case the wave 
function which solves the corresponding NLSE acquires the 
form 

^(x, t) = 2i/^e^^{x[l + Aexp(-^2)](1 + C')}"'/', (24) 
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FIG. 6: (Color on line) Plots of |^(x,t)|2 with A = 0.5 for the 
periodic (a) or quasiperiodic (b) bright soliton, with a = 0.1 and 
(3 = or a = (3 = 0.1, respectively, in the case of cubic and quintic 
nonlinearities. 
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FIG. 7: (Color on line) Plots of |^(x,t)|2 with A = 0.5 for the 
periodic (a) or quasiperiodic (b) dark soliton, with a = 0.1 and P = 
Oora = (3 = 0.1, respectively, in the case of cubic and quintic 
nonlinearities. 

where (j) = (j){x ^ t) is real, obtained via the Eq. (TTbl) . 

To get to the darklike solution we choose, for instance, fi = 
3, Gs = 6, and G5 = — 3. In this case the potential becomes 

^(x,t)=c^(t)x2 + /(x,t)-(3/4)l/\-'[l + Aexp(-e2)]^ 

(25) 

with ^3(x, t) and ^5(x, t) given by ^ and (l20l) , respectively. 
Moreover, for /i = 3, and for G3 = 6 and G5 = —3, we can 
use Eq. ^ to get to the solution ^{Q = C/\/l + It has 
the form of a kink, and so |<^>(C) P is now of the dark type. In 
this case, the solution given by Eq. (|5]) takes the form 

=6i/V*C{x[l + Aexp(-e2)](l + C')}-i/2. (26) 

We now illustrate the presence of bright and dark solitons in 
the case of cubic and quintic nonlinearities. In Figs. 6(a) and 
6(b) we depict according to (l24l) . for u;{t) periodic and 
quasiperiodic, respectively. In Figs. 7(a) and 7(b) we depict 
I^P according to (l26l) , for uj{t) periodic and quasiperiodic, 
respectively. Here we also note that the solutions evolve in 
time periodically or quasiperiodically, depending on the way 
we choose x{t) in (fT2l) , with no further changes, showing that 
they are stable localized excitations which we name periodic 
and quasiperiodic bright and dark solitons for the correspond- 
ing cubic and quintic NLSE. 

In summary, in this work we have obtained explicit soliton 
solutions of both the bright and dark type in the case of cubic 
and cubic plus quintic nonlinearities, in the presence of poten- 



tials which oscillate in time from attractive to expulsive behav- 
ior. The procedure is valid for a diversity of nonlinearities, and 
the results lead to the case of periodic and quasiperiodic so- 
lutions, having bright and dark features which evolve in time 
maintaining their localized behavior. As we have shown ex- 
plicitly, the system can support bright and dark solitons in the 
cases of cubic, and cubic plus quintic nonlinearities. The free- 
dom to choose the width of the solution has allowed to de- 
scribe periodic and quasiperiodic solutions, which are of cur- 
rent interest to the study of nonlinear excitations in conden- 
sates. Other possibilities of practical interest follow naturally. 

The authors would like to thank CAPES, CNPQ and 
PRONEX-CNPq-FAPESQ for partial financial support. 
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